Dissipation-driven two-mode mechanical squeezed states in optomechanical systems 
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In this paper, we propose two quantum optomeclianical arrangements tiiat permit the dissipation- 
enabled generation of steady two-mode mechanical squeezed states. In the first setup, the mechanical 
oscillators are placed in a two-mode optical resonator while in the second setup the mechanical oscil- 
lators are located in two coupled single-mode cavities. We show analytically that for an appropriate 
choice of the pump parameters, the two mechanical oscillators can be driven by cavity dissipation 
into a stationary two-mode squeezed vacuum, provided that mechanical damping is negligible. The 
effect of thermal fluctuations is also investigated in detail and show that ground state pre-cooling 
of the oscillators in not necessary for the two-mode squeezing. These proposals can be realized in a 
number of optomechanical systems with current state-of-the-art experimental techniques. 

PACS numbers: 



I. INTRODUCTION 

Quantum squeezing and entanglement have been ob- 
served in a number of atomic and photonic systems and 
are expected to play an increasing role in applications 
ranging from measurements of feeble forces and fields 
to quantum information science [H-Q- For example, it 
has been know for over three decades that squeezed vi- 
brational states are of importance for the measurement 
beyond the standard quantum limit of the weak signals 
expected to be produced e.g. in gravitational wave an- 
tennas dj. 

Although achieving such effects in macroscopic sys- 
tems remains a major challenge due in particular to the 
increasing rate of environment-induced decoherence @, 
recent progress toward th e g round-state cooling of mi- 
cromechanical systems |7Hllj may change the situation 
significantly in the near future. In particular, the char- 
acterization of quantum ground-state mechanical motion 
[l^ , the quantum control of a mechanical resonator deep 
in the quantum regime and coupled to a quantum qubit 
[l3| , and the demonstration of optomechanical pondero- 
motive squeezing are important steps toward the 
broad exploration of quantum effects in truly macro- 
scopic systems [T5l - [23| . Several proposals have been put 
forward to generate mechanical squeezing in an optome- 
chanical oscillator, including the injection of nonclassi- 
cal light [IJ], conditional quantum measurements [25| . 
and parametric amplification [26l - [29j . In all cases, deco- 
herence and losses are a dominant limiting factor in the 
amount of squeezing that can be achieved. 

A new paradigm in quantum state preparation and 
control has recently received increased attention. Its key 
aspect is that it exploits quantum dissipation in the gen- 
eration of specific quantum states. Quantum reservoir 
engineering has been proposed to prepare desirable quan- 
tum states [sol - lssj and perform quantum operations [s^ , 
and the creation of steady-state entanglement between 
two atomic ensembles by quantum reservoir engineering 
has been experimentally demonstrated [s^ . This dissipa- 
tive approach to quantum state preparation presents the 



double advantage of being independent of specific initial 
states and of leading to steady states robust to decoher- 
ence. 

Here, we propose to exploit cavity dissipation to gener- 
ate the steady-state two-mode squeezing of two spatially 
separated mechanical oscillators. These states are also 
entangled states of continuous variables, a basic resource 
in quantum information precessing. We propose specifi- 
cally two different setups. In the first one two mechanical 
oscillators are placed inside a two-mode optical resonator, 
while in the second one the mechanical oscillators are lo- 
cated in two separate single-mode cavities coupled by 
photon tunneling. The cavities are driven in both cases 
by amplitude-modulated lasers. In the first model we 
show analytically that for appropriate mechanical oscil- 
lator positions and pump laser parameters the mechani- 
cal oscillators can be driven into a stationary two-mode 
squeezed vacuum by cavity dissipation, provided that me- 
chanical damping is negligible. In the second case a two- 
step driving sequence can likewise give rise to a steady 
two-mode mechanical squeezed vacuum state. The effect 
of thermal fluctuations on the resulting squeezed states 
is also investigated in detail. 

The remainder of this paper is organized as follows. 
Section II introduces our first model of two mechanical 
oscillators coupled by radiation pressure to the field of a 
common two-mode optical resonator. It then discusses 
how to prepare a two-mode mechanical squeezed state. 
Section III shows that cavity damping can likewise be 
exploited to create a squeezed mechanical state of two 
mechanical oscillators in separated single-mode cavities. 
Finally, Section IV is a conclusion and outlook. 



II. MECHANICAL OSCILLATORS IN A 
SINGLE TWO-MODE CAVITY 

A. Model and equations 

We first consider an extension of the "membrane-in- 
the-middle" arrangement of cavity optomechanics (36j 



where two mechanical oscillators, modeled as vibrating 
dielectric membranes of identical frequencies are lo- 
cated inside a driven two-mode optical resonator, as il- 
lustrated in Fig. 1. The mechanical modes are charac- 
terized by the bosonic annihilation operators Cj and the 
cavity modes, at frequencies uJcj (j = 1,2), by annihila- 
tion operators Aj . These modes arc driven by amplitude- 
modulated lasers of frequencies ujij. The Hamiltonian of 
the driven cavity-oscillators system reads {h = \) 

H = ^ [ ^^cj a] Aj + LUm C] Cj + gjkA] Aj (Ck+Cl) 

+ i£,{t)e-'^'^'A]~t£*{t)e:^'^'AjY (1) 

where £j{t) are the time-dependent amplitudes of the 
pump lasers. Their specific forms will be given later. 
The single-photon optomechanical coupling constants 
= "X'^m^T""' HI' where L is the length of the cav- 
ity, m are the identical masses of the membranes, and 
fjkixk) = /f^^a'^^^i'A'r''- ^ ^ Here Uk and Xk arc the 
reflection coefficients and equilibrium positions of the 
two membranes and fccj the wave numbers of the cav- 
ity modes. 

For an appropriate combination of cavity length and 
membrane positions of the membranes, see Fig.l, it is 
possible to find a situation such that the "symmetrical" 
and "antisymmetric" optomechanical coupling strengths 
satisfy the equalities 

511 = .912 = .91, 

521 ^ -522 = .92, (2) 

as discussed in Ref. [s^l- Introducing then the new 
bosonic annihilation operators 

Bi = (Ci + C2)/V2, (3a) 
B2 = (Ci - C2)/V2, (3b) 

the Hamiltonian H becomes H = J^j Hj , where 

H, = Lo,,A]A, + u,r.B]B, + g.jA]A,{B, + b]) 

+ iEj{t)e~''^'^'A]~iE*{t)e"^'^*Aj. (4) 

In terms of the normal operators Bj , the system is there- 
fore decoupled and reduces to two independent single- 
membrane optomechanical systems. 

Further decomposing the operators Aj and Bj as 

= {A,{t))+a,=aj{t) + aj, 
B, = {B,it))+b,^p,it) + b,, (5) 

and with |aj(t)p > {a]aj) and |/3j(f)p > {b]bj) one can 
linearize the Hamiltonians Hj to get 

i/j'" - A,a]a,+c.,J]fe, + [x*(t)%+X,(Oa]](fc,+fcl), (6) 



2 




FIG. 1: Schematic plot of two vibrating membranes (Cj) 
placed at appropriately-chosen positions in a driven cavity 
with two frequency- nondegenerate resonant modes (Aj). 

where 

A, =<5, +.g, (/?,+/?*), 

with 5j = LOcj — coij, and the effective optomechanical 
coupling strengths are given by 

Xj{t) = gjUjit). (7) 

The density matrix pj of the sub-system composed of 
the cavity mode aj and the normal mode bj satisfies the 
master equation 

Pj{t) = -i[Hl''',p.j]+'^{2ajPja]~a]ajPj-pja]aj) 

+ ^(^th + mb^pfb] - bib.p, - p/bib,) 

+ -^"th(2&]pj&j - bjb]pj - pjbjb]), (8) 

where Kj is the cavity dissipation rate and 7„i the me- 
chanical damping rate, taken to be the same for both 
oscillators. The mean thermal phonon number is fith = 
{jJi^rn/kBT _ with Ub the Boltzmann constant and 

T the temperature. 



B. Stationary two-mode mechanical squeezed 
vacuum via cavity dissipation 

We now show how cavity dissipation can be exploited 
to prepare the stationary two-mode mechanical squeezed 
vacuum of the vibrating membranes. To this end we 
consider an effective optomechanical coupling strength 
Xj (t) of the form 

xAt)-Xjie-'^''^''^^^+X,2, (9) 

where Xji ^'^d Xj2 sre constants. This situation can be 
realized by a pump laser of the form discussed in Sec- 
tion ITCl 

For weak optomechanical coupling we have Aj ~ Sj , so 
that if the cavity-laser detuning Sj and the modulation 
frequency flj are 

Sj = Um, = 20;™, (10) 
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the transformation aj a^e and bj bje 
reduces the Hamiltonian to 



H 



lin 



(11) 



This form can be further simpHfied for a mechanical fre- 
quency iJrn S> Xjk, in which case the rapid oscillating 
terms e*^*"'"* can be neglected and 



{Xjie-'^^b,+Xj2b])aj+R.c. 



(12) 



This Hamiltonian describes the coupling of the cavity 
field Uj to the normal mode bj simultaneously via para- 
metric amplification - through the term proportional to 
Xji~ and via a beam splitter-type coupling - through the 
constant contribution Xj2 ■ It is known that the first term 
in the Hamiltonian ([T2|) leads to photon-phonon entan- 
glement and phononic heating of the normal mode bj, 
while the second term results in quantum state transfer 
between the cavity mode aj and the mechanical mode 
bj as well as to mechanical cooling (cold damping). To 
ensure the stability of the system, the coupling strengths 
must satisfy the inequality Xj2 > Xjii that is, cooling 
should dominate over anti-damping. 

When absorbing a laser photon the parametric am- 
plification term results in the simultaneous emission of 
a photon into the cavity mode aj and a phonon in the 
normal mode bj , while the beam-splitter interaction cor- 
responds to the annihilation of a photon and the emission 
of a phonon. We now show that the combined effect of 
these processes is the generation of steady-state single- 
mode squeezing of the normal mechanical mode bj, or 
equivalently of two-mode squeezing of the original modes 
ci and C2. 

We proceed by first performing the unitary transfor- 
mation 



Q3{t) = m,yp,{t)s,{^,), 



where the squeezing operator 



5(0) -exp -e*&f/2 + e,6|/2 



with 







(13) 



(14) 



(15) 

= tanh"^(xji/xi2)- (16) 
For 7,„ = the master equation ([5]) becomes then 



Qj{t) = -i\Hf\ Qj] + '^{2ajQja\ - a\ajQj - Qja]aj), 

(17) 

(18) 



where 



= G,(a]6j+6]a,), 



with 



Gj = X32\ll~(xrJxj^. 



(19) 



describes quantum-state transfer between the cavity 
mode Oj and the mechanical mode bj in the transformed 
picture. It can be inferred from that master equation 
that in the transformed picture and for 7^ = the cav- 
ity mode aj and the mechanical mode bj asymptotically 
decay to the ground state 



e,(00) = |0a,0fc^)(0a,0b^. 



(20) 



Simply reversing the unitary transformation (|13p we have 
then that in the steady-state regime the normal mode bj 
is indeed in the squeezed vacuum state 



P6^.(oo) = 5(0)|0,,)(0fc,|^n0)- 



(21) 



It is possible to adjust the amplitude Xjk and phase (pj 
of the optomechanical coupling coefficient ([9]) in such a 
way that the squeezing parameters satisfy the conditions 



r = Tj, 

01 = 02 - TT, 



(22a) 
(22b) 



in which case the two normal modes bi and 62 exhibit 
the same amount of steady-state squeezing, but in per- 
pendicular directions. With Eqs. ^ and (PT|) wc then 
have 

PccA^) = 5i2(62)|0e,,0c,)(0c,,0eJ^l2(62) (23) 

where we have introduced the two-mode squeezing oper- 
ator 



S'i2(62) = exp(-^*2Ci4 +62C1C2) 



(24) 



and £,12 ~ re^*"^. This two- mode squeezed vacuum of 
the mechanical oscillators can be thought of as the out- 
put from a 50:50 beam splitter characterized by the uni- 
tary transformation ([3]), the two inputs being the nor- 
mal modes bj squeezed in perpendicular directions. This 
shows that for 7^ — > the dissipation of the intracavity 
field can be exploited to prepare pure two-mode mechan- 
ical squeezed vacuum state. The minimum time imin re- 
quired for preparing such states can be evaluated from 
the eigenvalues of Eq pT|) , 




For symmetric parameters k 



one finds in 



4/k for G > k/2. 



Xj = Xjk, and Gj = G 



C. The choice of pump lasers 



Next we consider the choice of pump laser amplitudes 
£j{t) that result in the time-dependent optomechanical 
coupling © , and the specific form ([^ that results in the 
pure two-mode mechanical squeezing (P^ . From Eq. ([7]) 
an obvious starting ansatz is 



(25) 
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where ipjk are the initial phases of two components with 
amplitudes f and Jlj are modulation frequencies. The 
corresponding amplitudes aj{t) and Pj{t) of the cavity 
and mechanical modes are 



di 



/3j W = -(7m + i^^m)f3j - iffjiajf • 



(26) 



It is difficult to find exact solutions of these equations in 
general. For the case of the weak optomechanical cou- 
pling strengths gj^ however, approximate analytical solu- 
tions can be found by expanding the amplitudes aj and 

(3j in powers of gj as aj = aj"'' + a'^^^ + a'"^^ + ■ ■ ■ and 
= /3]°' + /3]i' + + . . . . Substituting these into 
Eqs. and for the time cj^ ^ 7m, = ^m., 

and flj = 2ujm one then finds (higher-order corrections 
can be derived straightforwardly) 



a 



(0) 

j 

(1) 



£72 



a,--' - 0, /3f 0, f3f ^ 0, 



„(2) 



2 ^2 



(27a) 
(27b) 



+ 



3w„i(k2 -I-w^J 
2zg|g,-i(3g^i+2£|,) 



3Wm(K^ + W2„)2 

2ig]£ji£j2 



-i[njt-{2<f,j-ipji)] 



3a;,„(K2 + w2j2 



„i(f2jt-ipji) 



2ig'^£l£j2 



Scomin'^ + a;^)^/^(K - 3iiu,„) 



-2i{Qjt-4>j) 



, (27c) 



^ a;™(K2+w2^) 3a;™(K2 -l-iu^) 



_l_ gj£ji£j2 ^ 
with the phases 



-i{Qjt — (j}j) 



't>3 = Vjl+V32. 

Lpj2 = arctan(a;„i/K). 



(27d) 



(28) 



For coupling strengths gj <C {k, e^fe, aj,„}, for example 
for gj ^ \Q~^ijJrm K ^ O.Q5oj„i, and £jk ^ lO^Wm, we find 



\af\ ~ 1 « ~ 10^ and S, ~ c^,„ » 5,(/3, +/?*)- 

10~'*a;m. We can therefore set Aj ~ Sj + gj{Pj + /3|) — 
5j, and the effective optomechanical coupling strength 

Xj{t) — 9jC<^f^ takes the required form © for 



JO) 



Xjk 



gj£jk 



(29) 



Finally, Eqs. ^ and (f22b| give 

^21 — 'fill — TT, ipj2 = arctan(wm/K). (30) 



D. Thermal fluctuations 



For finite mechanical damping, 7,„ ^ 0, the mechanical 
modes ci and C2 are no longer be in a pure squeezed state, 
but rather in a two-mode squeezed thermal state. In the 
state, we find from Eqs. ^ and ^ 



\CjCj) 00 
{ciC2)oa 



dodi cosh 2r — dod2 sinh 2r + sinh^ r, (31a) 

[- {dodi + 1 /2) sinh 2r + 4^2 cosh 2r] e"^, 

(31b) 



where di — nth cosh2r + sinh^ r, (i2 = (rith + 5) sinh2r, 
and 



rfo = 1 - 



('« + 7m)(K7m +4G2)' 



(32) 



The quantum correlations between the mechanical os- 
cillators can be quantified by the sum of variances [3, 113 



A 



EPR 



= ((xf^ +1,^^)2) + -r,^+^r), (33) 



where Xj^ = {cje~^^^ -|-cje*^^)/-\/2 arc quadrature opera- 
tors with local phase 6j. A value of Aepr < 2 is a signa- 
ture of Einstein- Podolsky-Rosen (EPR)-type correlations 
between the two mechanical modes, with Aepr = cor- 
responding to the ideal quantum mechanical limit 0] . In 
our case we find that AEPR.min is minimum for the choice 
of local phases Oi + 62 = 4>- In the long-time limit it is 
equal to 



iEPR,min 



2e-^''(l - do) + 2(2nth + l)do. 



(34) 



We first observe that for k = we have do = 1 and 
hence Aepr (00) = 4r;,th + 2, showing that in the absence 
of cavity dissipation the oscillators arc just prepared in 
the thermal state imposed by their mechanical coupling 
to the environment. This confirms that cavity dissipation 
is an essential component of this realization of steady- 
state two-mode mechanical squeezing. From Eq. p4p . 
steady-state squeezing requires that 



n-th < '^th.max 



1-rfo 
2do 



1 



(35) 



In practice, it is important to be able to operate at as high 
a number of mean thermal phonons as possible. This can 
be achieved by decreasing do while keeping the "cooper- 



ative parameter" 



4G^ 

K7„ 



» 1. From Eq. and for the 



realistic case k ^ 7™, do reduces approximately to 



do 



~ 4G2 ' 



(36) 



indicating that by increasing the coupling frequency G 
while keeping the ratio X1/X2 fixed it is possible to in- 
crease the value of fith max which is approximately given 
by 



nth,: 



4kxi(x2 - Xl) 
7™ [k2 + 4(x1 - X?)] ■ 



(37) 
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FIG. 2: (a) Two vibrating membranes (Cj) in separate single- 
mode cavities (^4^) are optically coupled by either an optical 
fiber, (b) Two evanescent-wave-coupled optomechanical crys- 
tal nanocavities. 



As a concrete example, for a cavity dissipation rate of 
K O.OSwm, a mechanical damping 7^ ^ 10~^a;m, and 
coupling strengths xi ~ 0.01a;„i and X2 ^ 0.03a;„i, we 
have fith,max ^ 70. This indicates that two-mode me- 
chanical squeezing is robust against thermal fluctuations 
and ground-state precooling of the mechanical modes 
may not be necessary. 



III. MECHANICAL OSCILLATORS IN 
SEPARATE SINGLE-MODE CAVITIES 

In this section we show that the same dissipative 
approach can also be utilized to generate a two-mode 
squeezed state of two mechanical oscillators in sepa- 
rate single-mode cavities. The specific example that 
we consider consists of two identical single-mode cavi- 
ties of frequency uJc, optically coupled to via fiber (a) 
or evanescent- wave coupling (b), see Fig. 4. Each cavity 
field is driven by a modulated laser and optomechani- 
cally coupled to a mechanical oscillator of frequency w,„ . 
Adopting the same symbols as before the Hamiltonian of 
that system reads 

3 =1,2 

+ Ji2{AiAl+A\A2), (38) 

where g account for the optomechanical coupling 
strengths, assumed to be identical, and J7i2 describes the 
coupling between the two cavities. Here we consider two 
pump fields of equal frequency and time-dependent 
amplitude 

£{t) = f ie"*"*+*'^i + £2e"^\ (39) 

For that symmetrical situation, both cavity fields have 
the same amplitude aj(t) = a{t) and the linearized 



Hamiltonian is 

Hiin = Y. + + + X{t)a]]{cj + c]) 

3 

+ Ji2{aial+a\a2), (40) 
where 

X{t) = ga{t) (41) 

and the detuning A ~ (lOc^'^i) for the weak optomechan- 
ical coupling. By introducing the new bosonic operators 

6i,2 = (ci±C2)/y2, (42a) 
di,2 = (ai ± a2)/V2, (42b) 

the Hamiltonian separates as before into the sum of two 
uncoupled Hamiltonians, Him = X]j=i 2 ^j™' where 

Hf- = A,d]d, + ioJlb, + [x*{t)d, + x{t)d]]{b, + S]), 

(43) 

with effective detunings Ai = A-|- J12 and A2 = A — J12. 
The dynamics of the two independent optomechanical 
sub-systems are governed by the same master equations 
as before, see Eq. ([5]). 

In contrast to the preceding case, though, in the two- 
cavity setup it is not possible to simultaneously achieve 
the single-mode squeezing of the normal modes bi and 62 
since there is only one driving laser field and one set of 
control parameters (wj, J7, (/3j). However, for sufficiently 
weak mechanical decoherence it is possible in principle to 
implement a two-step process that can still achieve that 
goal. 

For the first step we choose the frequencies lui and 57 
such that 

Ai = — + i/l2 = Wm.) 

n = 2l^™, (44) 

and the phases 

ipi = 01 — arctan(a;„j/K) 

ip2 = arctan(a;„j/K), (45) 

where 01 is arbitrary. With the transformation aj — >■ 
a^e^^^^* and bj — > 6^6"'""*, the Hamiltonians ifj'" re- 
duce to 

H'r = {xie-"^'h+X2b\)di 

+ (X2e~'''"'"*6i + xie'(2"'"*-'^^^&l)di + H.c. 

(46a) 

i?2" = (xie^^'-^"*-^^^^i +X2e-^'-^-*6l)Ji 

+ (xie''^-^"-"")*^: + X2e''[2(-^-+'^'"^*-^'^^l&I)di 
+ H.c. (46b) 
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For Xj ^ {^m, J^i2, \ J\2 — the non-resonant terms 

in these Hamiltonians can be neglected and they reduce 
to H}'" ~ (xie^^'^^Si + X2h\)k + H.c. and Hf" ~ 0, 
respectively. For the mode h\ this is formally the same 
situation as encountered in Section IIB. Neglecting as 
before the mechanical damping, 7„i = 0, cavity dissipa- 
tion brings likewise that normal mode into a steady-state 
single-mode squeezed vacuum for long enough time, and 
at the same time mode hi simply decays into the vacuum, 
i.e., 

~PlS^ > imin) = ^(Ci)|0,,)(0,J5+(a), (47a) 

/5^,(i>imin)- |06,)(0bJ. (47b) 

where imin, as defined before, is a time long enough that 
steady state has been reached. After the first step the 
mechanical oscillators ci and C2 are therefore prepared 
in a pure two-mode squeezed state, although not a stan- 
dard two-mode squeezed vacuum. That latter goal can 
be achieved in a second step by changing the frequency 
and phase of the pump laser at time t„ii,i so that 



with 



A2 

n 



UJc - UJl 



Jr. 



and the phases 

V2 



(f>i — arctan(a;„i/K) + tt, 
arctan(a;„i/ n). 



(48) 



(49) 



0th,bi 
f 



(r — f)e^"^, 

1 2do(di+rf2) + l 

4 °2do(di-d2) + l' 



— |nbj(nbj, i 



1,2 



(53) 



with the mean thermal excitation number ni = nth and 
fi2 = y^{dodi + 1/2)2 _ ^2^2 _ 1^2. Note however that 
a two-step procedure is not required in that situation 
since a single step with laser parameters satisfying either 
Eq. [ pS)) ] will result in the preparation of mode bi 
[or bi] in a squeezed thermal state while the other mode 
in the thermal state (j52ap . with the degree of EPR cor- 
relations between the mechanical modes ci and C2 



A 



EPR, mill 



(oo) = e-2-(l - do) + (2nth + 1)(1 + do). 



(54) 



Hence, the maximum mean number of thermal phonons 
TT-th max for which squeezing can be achieved is 



1 



^th,max 



2(1 + do) 



(1-e-n 



(55) 



which is obviously smaller than that in Ea. p5|) because 
just one normal mode is now in a squeezed state and it 
is therefore harder to maintain quantum correlations. 



For the weak optomechanical coupling, we then have that 
jjhn ^ which means that the mode bi evolves freely 
in the state of Eq. (|Tfal) . while H^"^ ~ (xie~*'^'&2 + 



X2bl)d2 



H.c. After a time t > 2t, 



both the normal 



modes bi and 62 are therefore in steady-state single-mode 
squeezed vacua, 

Pb^{t > 2t,„i„) = ^(6)|06j(06j5t(ei), (50a) 

Pb,{t > 2tmi„) = 5(6)|06,)(0bJ^t(^2), (50b) 

and the mechanical modes ci and ci are in the two-mode 
squeezed vacuum state. 



3 (^ ^ 2tj^in) 



^12(a2)|0e,,0c.)(0ci,0c,|5l2(a2). 



(51) 

When accounting for mechanical damping the two- 
step preparation scheme remains efficient for mean ther- 
mal phonon number such that 7m ^th ^ k so that 
[7m^th]~^ ^ imin ^ and thermal effects can be 
neglected during the state preparation. In case this con- 
dition is not satisfied, 7m"-th ^ i^, following the second 
step the mode bi is thermalized while the mode 62 is in 
a squeezed thermal state, 

/Obi(t > 2t,„j„) = pth,&i, (52a) 

PbAt > 2im^n) = Si^2)pth,bJH^2), (52b) 



IV. CONCLUSION 

In conclusion, we have proposed two possible quantum 
optomechanical setups to generate two-mode mechani- 
cal squeezed states of mechanical oscillators in optical 
cavities driven by modulated lasers. We showed analyti- 
cally that for appropriate laser pump parameters the two 
oscillators can be prepared into a stationary two-mode 
mechanical squeezed vacuum with the aid of the cavity 
dissipation when choosing appropriate positions of the 
membranes in the optical cavities. The effect of ther- 
mal fluctuations on the two-mode mechanical squeezing 
was also investigated in detail, and we showed that me- 
chanical squeezing is achievable without pre-cooling the 
mechanical oscillators to their quantum ground states. 
The present schemes are deterministic and can be im- 
plemented in a variety of optomechanical systems with 
current state-of-the-art experimental techniques. 
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